RHOMBUS TILINGS OF A HEXAGON WITH TWO TRIANGLES 
MISSING ON THE SYMMETRY AXIS 



THERESIA EISENKOLBL 



Abstract. We compute the number of rhombus tilings of a hexagon with sides n, n, N, n, 
n, N, where two triangles on the symmetry axis touching in one vertex are removed. The 
case of the common vertex being the center of the hexagon solves a problem posed by Propp. 



1. Introduction 

The interest in rhombus tilings has emerged from the enumeration of plane partitions in a 
given box. The connection comes from looking at the stacks of cubes of a plane partition from 
the right angle and projecting the picture to the plane. Then the box becomes a hexagon, 
where opposite sides are equal, and the cubes become a rhombus tiling of the hexagon where 
the rhombi consist of two equilateral triangles (cf. 0). The number of plane partitions in a 
given box was first computed by MacMahon |7|, Sec. 429, q — > 1, proof in Sec. 494]. Therefore: 

The number of all rhombus tilings of a hexagon with sides a, b, c, a, b, c equals 

i=lj=l fc=l J 

(The form of the expression is due to Macdonald.) 

In ||], Propp proposed several problems regarding "incomplete" hexagons, i.e., hexagons, 
where certain triangles are missing. In particular, Problem 4 of jg] asks for a formula for the 
number of rhombus tilings of a regular hexagon, where two of the six central triangles are 
missing. We treat the case of the two triangles lying on the symmetry axis and touching in 
one vertex (see Figure ||) . The other case has been solved in || . We prove the following two 
theorems. 

Theorem 1. The number of rhombus tilings of a hexagon with sides n,n,2m,n,n,2m and 
two missing triangles on the horizontal symmetry axis sharing the (s + l)-th vertex on the 
axis (see Figure \^) equals 

/o™ i A (2m— 2\ (2n— 2s\ /2s\ n n 2m . . , , 
(2m — 1) ( m _ 1 ) ( n _ s ){ s ) -i-r -i-r t-t I + J + k — 1 

/2m+2n\ 11 11 11 i , . k _ 2 

V m+n J i=lj=lk=l J 

Theorem 2. The number of rhombus tilings of a hexagon with sides n, n, 2m + 1, n, n, 2m + 1 
and two missing triangles on the symmetry axis sharing the s-th vertex on the axis equals 

(2m + l)( m )( n _ s )( 8 _J * + J+fc~l 

/2m+2nN nn n i+j + k _ 2 - 

\ m+n J i=\j=\ k=l J 

The following corollary is easily derived using Stirling's approximation formula. 

Corollary. The proportion of rhombus tilings of a hexagon with sides at, at, f3t, at, at, (5t and 
two missing triangles on the horizontal symmetry axis touching the (jt)-th vertex on the axis 
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Figure 1. A hexagon with sides n,n,2m,n,n,2m and missing triangles in 
position s + 1, where m = 2, n = 3, s = 2, and a rhombus tiling. 



in the number of all rhombus tilings of the hexagon with sidelengths at, at, fit, at, at, fit (given 
by B (at, at, fit, at, at, fit) in ffity) is asymptotically equal to 

J_ l fi(2a + fi) 
4-7T y 7(0 — 7) 

This expression can attain arbitrary large values if a is close to 7. It equals (which is 
approximately 0.28) for a = fi = 2j, which corresponds to the case of a regular hexagon with 
two missing triangles touching the center. In comparison, in the other case of Problem 4 of 
H , the case of a fixed rhombus on the symmetry axis, the analogous proportion must clearly 
be smaller than 1 and equals approximately | if the central rhombus is missing (see ||). 

The rest of the paper is devoted to the proof of Theorems [l] and |2[ The main ingredients 
are the matchings factorization theorem by M. Ciucu jjj], nonintersecting lattice paths, and 
two determinant evaluations, the latter constituting the most difficult part of the proof. An 
outline of the proof is given in the next section. The details are filled in in the subsequent 
sections. 



2. Outline of the proofs of Theorems [l] and || 
Outline of the proof of Theorem |]: 

Step 1: It suffices to compute the number of perfect matchings of two hexagonal graphs G + 
and G~ (see Section 

We use the fact that every rhombus tiling corresponds to a perfect matching of the (inner) 



dual graph (see Figure 
theorem by M. Ciucu 



a) . This graph has reflective symmetry, so the matchings factorization 
IJ is applicable (see Lemma |3|) . This theorem expresses the number of 



perfect matchings of a graph as product of 2™ and the numbers of perfect matchings of two 
smaller graphs G + and G~ (see Lemma |I| and Figure ^) , which are roughly the two halves of 
the original graph. The remaining task is to count the numbers of perfect matchings of G + 
and G~ . 

Step 2: The numbers of perfect matchings of G + and G~ equal the numbers of rhombus 
tilings of two regions R + and R~ consisting of triangles, respectively (see Sections and 
We convert the graphs G + and G~ back to regions R + and R~ consisting of triangles (see 
Figure |3|), so that we have to count rhombus tilings again. 

Step 3: The numbers of rhombus tilings of R + and R~ are certain determinants (see 
Sections [| and[^J. 

The rhombus tilings are in bijection with certain families of nonintersecting lattice paths 
(see Figures ^| and ||) . Application of the main result of nonintersecting lattice paths expresses 
the desired numbers as determinants (see Lemmas [5] and ||). 

Step 4-' Evaluation of the determinant corresponding to G + (see Section^). 
The determinant corresponding to G + is evaluated using a lemma by Krattenthaler (see 
Lemmas |6| and ^) . 

Step 5: Evaluation of the determinant corresponding to G~ . 

We pull factors out of the determinant, so that we obtain a determinant whose entries 
are polynomials in m (see the proof of Lemma ||). Then, in Lemma ||, this determinant is 
evaluated by using the "identification of factors" method, as explained in @, Sec. 2]. The 
corresponding details are the subject of Sections || - ||. 

Step 6: A combination of the previous steps proves Theorem 

We substitute the results of Lemmas 0, || and |9] in Lemma ||] and obtain the following 
expression for the number of rhombus tilings of our original hexagon, 

2 ^_ 1 (H(n)) 2 (2n-2 S -l)!!(2 S -l)!! 
H(2n)(n - s)\s\ 

n—2 / \ min(fc,n-l-fc) n 

X ]"[ (2m + 2j-i)[J (m + k + -\ JJ (m + jfc)min(A ! +l,n-fc+l) i 

2<i<j<n k=l ^ ' k=0 

Here, H(n) stands for nr=o This can easily be transformed to the expression in Theorem [j], 
so the proof of Theorem [l] is complete. 

Outline of the proof of Theorem |2] (see Section |9|): 

Step 1' : It is enough to count the rhombus tilings of the two regions R + and R~ shown in 
Figure 

By the analogue of Step 1, we have to count the numbers of perfect matchings of two graphs 
G + and G~ . These numbers equal the numbers of rhombus tilings of the regions R + and R~ 
of triangles shown in Figure ^ (compare Step 2 in the previous proof). 

Step 2' : The number of rhombus tilings of R + . 

The number of rhombus tilings of R + (n, m) equals the number of rhombus tilings of R + (n — 
1, m), so it can be found using Lemma (see equation (|T^)). 

Step 3' .' The number of rhombus tilings of R~ . 

The number of rhombus tilings of R~(n,m, s) can be expressed as a determinant analo- 
gously to Step 2 in the previous proof. The determinant equals the determinant counting 
the number of rhombus tilings of R~(n — 1, m + 1, s — 1), so Lemma ^ can be used for the 
evaluation (see equation fll3])). 
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Step 4': A combination of the previous steps proves Theorem ||. We substitute the results 
of equations © and (||) in (0). We get 

M(G) = 2 n ~ 1 M(E+)M(R 1 ) 

_ 2( n 2 1 )" 1 H(n - 1) H(n + l)(2n - 2s - l)!!(2s - 3)!! 
(n _ s)! ( s _l )!n n =o(2j - )!n -2 (2z + 1)! 

n—2 , 1 v min(fc— l,n— fc-1) « 

x]^[fm + A: + -j JJ (m + jfc) min (*- n - fc+1 ) \\ (2m + 2j-i), 

k=2 ^ ' k=l 2<i<j<n+l 

which can easily be transformed to the expression in Theorem ||. 

3. Breaking the hexagon in two parts 

We start the proof of Theorem |l| by forming the inner dual of the given hexagon. I.e., we 
replace every triangle by a vertex and connect vertices corresponding to adjacent triangles. 
Thus, we get a hexagonal graph, whose perfect matchings correspond to rhombus tilings of 
the original hexagon (see Figure |2|a) . 

Now we use a theorem by M. Ciucu (see [||) to, roughly speaking, break the hexagonal 
graph into two halves. This theorem is described as follows. Let G be a graph with reflective 
symmetry, which splits into two parts after removal of the vertices of the symmetry axis. 
Label the vertices on the symmetry axis oi, b\, 02, 62, • • • , Oqi, 021 from left to right. If G is 
bipartite, we can colour the vertices of the graph black and white subject to the conditions 
that a% is white and no two adjacent vertices are of the same colour. 

Then we delete all edges connecting white a-vertices and black 6-vertices to the upper half 
and all edges connecting black a-vertices and white 6-vertices to the lower half. If we divide 
by two all weights of edges lying on the symmetry axis, the graph splits into two parts G + 
and G~ . Now we can state the matchings factorization theorem from [|J. 

Lemma 3. Let G be a planar bipartite weighted, symmetric graph, which splits into two parts 
after removal of the vertices of the symmetry axis. Then 

M{G) = 2 l ( G) M{G + )M(G~), 

where M(H) denotes the weighted count of perfect matchings of the graph H and G^ denote 
the upper and lower half of G as described above. 21(G) is the number of vertices on the 
symmetry axis. 

We apply Lemma ||| to our hexagonal graph, exemplified in Figure |2ja, with respect to the 
horizontal symmetry axis. In our case 1(G) = n — 1, since two vertices correspond to the 
removed triangles. G + and G~ are shown in Figure |2|b. Thus we get the following lemma. 

Lemma 4. The number of rhombus tilings of a hexagon with sides n, n, 2m, n, n, 2m and two 
missing triangles on the symmetry axis sharing the (s + l)-th vertex on the axis equals 

2 n - 1 M(G + )M(G~), 

where G + and G~ are formed by the above procedure, as exemplified in Figure [|&. 
M(G + ) and M(G~) are computed in the following sections. 

4. The matchings count for the upper half 

In this section we evaluate M(G + ). We start by expressing M(G + ) as the following deter- 
minant. 
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a. The perfect matching corresponding b. G + and G" 

to the rhombus tiling of Figure [|. 

Figure 2. 

Lemma 5. 

M(G+)= det (Y^" 1 

l<i,j<n \ \m — J +1 

Proof. First, we convert G + back to the corresponding region of triangles, R + say (see Fig- 
ure ||a), so that perfect matchings of G + correspond bijectively to the rhombus tilings of R + . 
Thus, we have to count rhombus tilings of R + . The next step is converting rhombus tilings 
to families of nonintersecting lattice paths, i.e. lattice paths which have no common vertices. 
The reader should consult Figure ||, while reading the following passage. Given a rhombus 
tiling of R + , the lattice paths start on the centers of upper left diagonal edges (lying on one of 
the sides of length n). They end on the lower right edges parallel to the starting edges. The 
paths are generated by connecting the center of the respective edge with the center of the edge 
lying opposite in the rhombus. This process is iterated using the new edge and the second 
rhombus it bounds. It terminates on the lower right boundary edges. It is obvious that paths 
starting at different points have no common vertices, i.e., are nonintersecting. Furthermore, 
an arbitrary family of nonintersecting paths from the set of the upper left edges to the set 
of the lower right edges lies completely inside R + and can be converted back to a tiling (see 
Figure |||a). 

Then we transform the picture to "orthogonal" paths with positive horizontal and nega- 
tive vertical steps of unit length (see Figure |]b,c). Let the starting points of the paths be 
denoted by Pi, P2, ■ ■ ■ ,P n and the end points by Q±, Q2, ■ ■ ■ , Q n - We can easily write down 
the coordinates of the starting points and the end points: 

Pi = (i — 1, i + m — 1) for i = 1, . . . , n, 

Qj = (2j-2,j-l) for j = 1, ... , n. 



Next we apply the main result for nonintersecting lattice paths [Q, Cor. 2] (see also [10, The- 
orem 1.2]). This theorem says that the weighted count of families of nonintersecting lattice 
paths, with path i running from P, to Qi, is the determinant of the matrix with (i, j)-entry 
the weight P(P« — > Qj) of lattice paths running from Pj to Qj, provided that every two paths 
Pi — > Qj and Pk — > Qi have a common vertex if i < j and k > I. It is easily checked that our 
sets of starting and end points meet the required conditions. 



a. R + , the upper half of the hexagon. 



b. R , the lower half of the hexagon. 



Figure 3. 

The number of lattice paths with positive horizontal and negative vertical steps from (a, b) 
to (c, d) equals ( c_ ^p d ) ■ Therefore, the number of families of nonintersecting lattice paths 
(equivalently, the number of rhombus tilings of R + ) is equal to the following determinant: 

det (V(P i ^Q j ))= det (( m + J ~ 1 
This proves Lemma [| □ 

This determinant can be evaluated with the help of the following determinant lemma ( || , 
Lemma 2.2). 

Lemma 6. 

det ((xj + a n )(xj + a n _i) • • • (xj + a i+1 )(xj + bi)(xj + . . . (xj + b 2 )) 

l<i,j<n 

= n ( x i~ x j) n ( b i~ a j)- 

l<*<i<« 2<i<j<n 

The determinant in Lemma [| factors as follows. 
Lemma 7. 



M(G + ) = M(R + ) 



H (™)n 2 < i < i <n(2m + 2j-i) 



ny=i(2j-2)i 

Here, in abuse of notation, M(R) denotes the number of rhombus tilings of R, if R is a 
region of triangles. 

Proof. Before we can apply Lemma ||, we have to transform the expression of Lemma |5| in the 
following way. 

+ J — 1 N 

m + j — 1)! 



M(G + ) = det 

l<i,j<n \ \m — j + i 



W 2)1 i n (n+m _ i)!(2i _ 2)! 



x det ( (m + i + 1 — j)(m + i + 2 — j) . . . (m + n — j) 

Ki,j<n \ 



-./ + -) ( + )•••(-./ + 1 



a. A tiling of the upper half of the hexagon 
and the corresponding lattice path family. 



b. The paths isolated. 



Pa 



Pi 



• ® 



Qi 



Q2 



Qs 



c. The corresponding lattice path family. 



Figure 4. 



Now we apply Lemma ^ with xj- 
result. 



-k, cik = m + k, bk = | and simplify to get the claimed 

□ 



5. The matchings count for the lower half 

In this section we reduce the computation of M{G~) to the evaluation of a determinant, 
which is stated in the following lemma. 



Lemma 8. 

M(G-) = M(R 



(n + m — s)(s + m) 



det(Bij), s = 0, . . . , n — 1 



(2n-2s)nr=i (2»+l-2t)! 
where R~ is the region of triangles which is described in the proof of the lemma, and where 

r (n + 2 + j - 2i) n _j(i + m + I - j)^ x (m + § + \-(\ i^s + l 
(n + l+j — 2s) n -j(s + m+ 1 — i = s + 1 



Bi 



(2) 



a. A tiling of the lower half of the hexagon 
and the corresponding lattice path family. 



b. The paths isolated. 



Ri 



R2 R3 

$ ®- 



S 3 



S 2 



Si 



c. The corresponding lattice path family. 



Figure 5. 



Here (a) n := a(a + l)(a + 2) 
% ^ s + 1 can also be written as 



(a + n — 1) is the usual shifted factorial. The entry for 



B 



1 



1.1 



-(n + 2i) n _ j+l (i +m + 1 - + (n + 2 + j - 2i) TV _ J -(i + m-j)j. 



(3) 



Proof. We start analogously to Section | and convert G~ (exemplified in Figure ^3) back 
to a region R~ of triangles (see Figure |b), so that the perfect matchings of G~ correspond 
bijectively to the rhombus tilings of Rr . However, since G~ contains edges on the symmetry 
axis of G, which, by Lemma [3l has the consequence that they count with weight | in G~ , we 
are dealing with a weighted count of the rhombus tilings of R~ , where rhombi such as the 
top-left rhombus in Figure ||a count with weight i. Again, we count the rhombus tilings of 
R~ by counting the number of nonintersecting lattice path families leading from upper left 
to lower right edges. The starting and end points can be easily read off Figure [5]. 



The starting points are (see Figure ||| for the meaning of s, 1 < s < n — 1): 

R;, = 



2, m + i — 1) for + 1 
• 1, m + s — 1) for i = s + 1. 



The end points are: 

£j = ( n + J - !>i - 1) for j = 1, . . . ,n. 

Now we apply again the main result for nonintersecting lattice paths. The matrix entries 
are V(Ri — > Sj). We note that a positive horizontal step starting at B4, i ^ s + 1 corresponds 
to a rhombus of weight \ in so paths starting with a horizontal step at R, i ^ s + 1 are 
counted with weight g. Therefore, we count paths starting with a horizontal step and paths 
starting with a vertical step separately. 



We get M(G ) = deti<jj< n (Ay), where 

Ay = -> Sj) 



' 1 / n + m — i\ ( n + m — i 
z \m + 1 — J J \m + 1 — 1 — j 



n + m — s 
m + s — j 



(4) 

for i = s + 1. 



Since this expression also makes sense for s = 0, we can include this case in the following 
calculations. We pull factors out of the rows of det(Ay), so that the remaining entries are 
polynomials in m and get the expression in Lemma ||. □ 

The next lemma states an evaluation of the determinant in Lemma J. An outline of the 
proof is provided immediately below. The details are given in Sections || - ^. 

Lemma 9. 

. . 2(V)H(ra)(2ra-2s-l)!!(2s-l)!! 

det (Bij) = — — — — — 

(n — s — l)!s!(m + s)(m + n — s) 



n—2 / . \ min(fe,n— 1— k) n 

x Y[(m + k + ^\ H(m + k) mi < k+1 ' n - k+1 \ 

k=l ^ ' k=0 



Outline of the proof of Lemma |]: In Section we prove that Ilfc=i { m + ^ + I) mm<yk,n 1 ^ 
divides det (Bij) as a polynomial in m. In Section ^ we prove that Ofc=o ( m fc) mm ( fc+1 ' n_fc+1 ) 
divides det (Bij) polynomial in m. 

The results are stated in Lemmas 10 and 11. We show this with the help of linear combi- 



nations of rows and columns, which vanish, if one of the linear factors is set equal to zero. 

In Section [|| we compute the degree of the determinant as a polynomial in m. It is exactly 
equal to the number of linear factors we have already found to divide the determinant. So we 
know the determinant up to a constant factor. We compute this constant in Lemma ^ by 
replacing each entry by its leading coefficient and using Vandermonde's determinant formula. 

Lemmas 10, 11 and ^ immediately give Lemma ||. □ 



6. The "half-integral" factors of det (Bij) 
In this section we prove the following lemma (see Lemma || for the definition of B^ ) . 

Lemma 10. Ilfe=i { m + ^ + l) mm ^ k ' n 1 k ^ divides det (Bij) as a polynomial in m. 



Proof. We find linear combinations of columns which give zero for m = — (k+ 2 )- First, we 
show that the following linear combination of columns equals zero for i^s + 1, 0<l<k, 
l> 2k-n + l: 



j=o 



D 



i,n+2l-2k-j 



0. 



(5) 



In order to establish this, we break the sum in two parts according to the two summands 
of Bij in equation (13) and convert them to hypergeometric form. The left-hand side of dBI) 



becomes 
(k-2l + i 



n) n _ 2k+2l (2-2k+2l-2i+2n) 2k _ 2l2 F 1 



1 



-1,-1 + 2k- 
2k-2l + i 



2l + 2i-2n 
l : 1 

2 



k 



n 



+ (l + 2k-2l + i-k--- n)_i_ 2 fe+2«+n(2 - 2k + 21 - 2i + 2n) 2k ^ 2i 



Now, Vandermonde's summation formula 



2^1 



-n 



2*1 



-a + c 



-l,-l + 2k-2l + 2i-2n 
l + 2k-2l + i-k- 



n 



(C)r 



is applicable to both 2-^i-series, since I > 0. It is directly verifiable that the two resulting 
expressions sum to zero. 

It is easily seen that the conditions on I in (S) allow min(/c + l,n — k) possible values for 
I. Thus, we have min(/c + l,n — k) independent linear combinations of columns which vanish 
in all coordinates except possibly in the (s + l)-th coordinate. (Recall that (|5|) is valid only 
for i 7^ s + 1). It is clear that an appropriate combination of two of these linear combinations 
vanishes in every coordinate. So we have min(fc, n — k — 1) independent linear combinations 
vanishing at m = — k — i for k = 1, . . . ,n — 2, which proves Lemma □ 



7. The "integral" factors of det (By) 
In this section we prove the following lemma (see Lemma || for the definition of By ) . 



Lemma 11. 



l 

(n+m— s)(s+m) 



YYk=o ( m + ^) min ( fc + 1 ' n k+l ) divides det (Bij) as a polynomial in 



m. 



Proof. We use linear combinations of the rows of By that vanish for m = —k. Without loss 
of generality, we can assume that s < l| because both the final result (see Theorem |l|) and 
the number of rhombus tilings of the original graph are invariant under the transformation 
s — > n — s. 

Most of the factors (m + k) can be pulled out directly from the rows of By. In fact, it is 
easily seen that row i is divisible by (m + 1 — i + n) 2 i-n-l f° r 2i > n + 2. The product of 
these terms equals Hfc=i ( m + /j) 111111 ^™-^) . The matrix (Cy) which remains after pulling out 
these factors from By looks as follows: 

(n + j + 1 — 2s) n -j(s + m + 1 — j)j-i i = s + 1, 

Cy = < (n + 2 + j - 2i) n -j(i + m + 1 — j)j-2i+n(2m + n + 1 - j) i / s + 1, 2i > n + 2, 
(n + 2 + j - 2i) n -j(i +m + 1 - j')j_i(2m + n + 1 - j) i / s + 1, 2i < n + 2. 

To finish the proof of Lemma 11 we have to find for each value k = 0, . . . ,n, k ^ s and 
k 7^ (n — s) one vanishing linear combination of the rows of Cy. 
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We start with the case k < s. We claim that 



E 

i=k+l 



-1) 



k - 1 



(n + 



-k-i(n + 



i-k-l 



- k - 1/ ( s + i - i) i _ k _ 1 (n - k-i + l)i-jfe_i 



2(n + | - s - l) s -k(n - s + 

(|) 5 -Jfe-l(n ~k-S + l) a _jfe_i 



'J 



G 



0. (6) 



If —A; + s — j < then the terms (i — k + 1 — .7)1—1 and (s — k + 1 — j)j—x (which are factors of 



and C s +i 



respectively) are zero for all occurring indices. If —A; + s — j > 0, 



we reverse the order of summation in the sum and write (||) in hyper geometric form, 



1) 



-k+s 



(n + 



s) s -k-i{l-j - k + s) 



3-1 



a) 



2 )s—k—l 



(n-k- s + l) s _ fc _i(2+ j + n -2s) n -j (n - s + l) s _k_i 



x (-1+ j + 2k-n) 3 F 2 



2,n — k + 1 — s,j + k — s 

1 + i + tt — S, j + i + if — ; 



;1 



+ (-i) 



s-fc+2 



2(n + | - s - l) s -k{n -s + l) s _fc-i 
(l)s-jfc-l(n - k - s + l) s _fe_i 



Gs+1, 



Now we can apply the Pfaff-Saalschiitz summation formula (II, (2.3.1.3); Appendix (III. 2)), 



3F2 



a, b, —n 
c,l + a + b — c — n' 



{-a + c) n (-b + c) n 
(c)„ (-a - 6 + c)„ 



(7) 



It is easily verified that the resulting sum of two terms equals zero. 
The case k > n — s is quite similar. We claim that 



E 

=n—k+l 



-1) 



i— n+k+1 



re + fc-l\(n+§- i)j_ n+fc _i(n - i + l);_ n+fc _i 



n + k-lj (s + i- i)j_ n+fc _i(A; - i + l);_ n+fe _i 



m=—k 



("I) 6 



-n+fc+2^( n + I — S — l)s-n+fc( n — S + l) s - n +fc-l 



(^)s-n+fc-l(^ — S + l) s - n +fc-l 

Converting the reversed sum to hyper geometric form gives 



a 



s+l,j 



0. (8) 



(-1) 



s+k—n 



(-1 +j + 2k - n)(2 + j + n- 2s) n -j (1 + n- s)_i +fc _ n+s 

3(1 — j — k + s)-i +j 



[% + n- s)_i+fc- 



l+fc— n+s 



(1 + fc - s)_i + jfc_ 



X 3^2 



1 — k + n — s, r;,j + k 



l+2~'~2 s i 2 ^ 2 2 S 



Again, the Pfaff-Saalschiitz summation formula (m) is applicable because — (1 — k + n — s) is 
a nonnegative integer. It is easily checked that the resulting terms sum to 0. 

So our remaining task is the case s < k < n — s. For s < k < S we consider the following 
linear combination, 
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E (- 4 )"" 

i=k+l 



(s-i + l)i_ fc _i(s - n + \ 



2 )n—k—l 



2n-2i + l)!(s + i - i)i_ fc _i(n + 1 - s)_ fc 

I n+l+j I 

- J - (s - i + l)j_ fe _i(s - n+ i 



+l-2iCij 



2 )n—k—l 



_L (-4)" v 

^3, (2n-2i + l)!( S + £-i) Wfc _ 1 (n + l-a)_ fc " 



-Ci 



Now, both the term (i — fc)n+l-2iCV 



8+1J 



0. (9) 



m=— fc 

which is part of the first sum in <M) and the term 



which is part of the second sum in (|9|) are equal to (n + 2 + j — 2i) n -j{i — k + 1 



j)j-2i+n(—2k + n + 1 — j), so we can combine the two sums into one sum of a hypergeometric 
term. 

We distinguish two cases according to the parity of n — j. In both cases we reverse the sum 
and convert it to hypergeometric form. The resulting two hypergeometric series are 



3 r 2 



3 F 2 



l + l + m,—l — m, ^ + I — n + s_ 



1 + l-n 



2' 

1 + / + m, — 1 — I — m, \ + I — n + s_ 



|, 1 + I — n + s 



for n — j = 21, 
for n — j = 21 + 1. 



So we can use the Pfaff-Saalschiitz summation formula @ again. It is easily verified that in 
both cases the resulting terms add to zero. 

The case ^ < k < n — s is handled similarly. We claim that the following sum equals zero, 



E (- 4 ) 

j=n— fc+1 

L 



(s-i + l)j_ n+fc _i(s -n + \)k-\ 



n+j + l 



+ E M 



(2n -2i + l)!(s + 5 - i)i- n+k -i{n + 1 - s)^ n+k 
{s-i + l)j_ n+fe _i(s - n + \)k-i 



(i - k) n+ i^2iQ 



m=—k 



(2n-2i + l)!(s + i-i) 



fc _i(n + 1 - s)_„ +fc 



m=— fc 



-l) n G 



m=— fc 



0. (10) 



Again, we write the two sums as one single sum, distinguish two cases according to the parity 
of n — j, and reverse the order of summation. Conversion to hypergeometric form of the 
resulting sums gives 



3-^2 



3-^2 



1 + I - k, -I + k, \ + 1 - n + s 
§,l + Z-n + s 

l + l-k,—l — l + k,^ + l — n + s 
|, 1 + I — n + s 



;i 



for n — j = 21, 
for n — j = 21 + 1. 



The Pfaff-Saalschiitz summation formula (^) can be applied in both cases. It is easily seen 
that the results vanish after subtraction of (— l) n C s +ij . Thus Lemma 11 is proved. □ 



m=—k 
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Figure 6. The hexagon in the case of odd sidelength and the two halves R + 
and R~ . m = 2, n = 3, s = 2. 

8. The degree and the leading coefficient 

We have to find the degree and the leading coefficient of the determinant det (Bij) as a 
polynomial in m. The degree of Bij is j — 1 for i = s + 1 and j else. The degree of det (By) 
is clearly ( n ^ ) — 1, which is easily seen to be the number of linear factors we have found 
to divide det (.Ay). Therefore det (By) is the product of the linear factors and the leading 
coefficient. 

To compute the leading coefficient we look at the leading coefficient of each entry. The 
leading coefficients of each entry give the matrix D, with 

, . (2-2i i^ s + 1 

Dij = (Xi + n + j) n -j, where x i = \i_ 2s i = s + l 

This matrix can be transformed by column reduction to (x™ -1 ), but this is just the Vander- 
monde determinant r[i<j<j<n ( x « ~ x j)- Plugging in the values of the xi gives the following 
lemma. 

Lemma 12. The leading coefficient of det (By) is 

2(V) H(n)(2n - 2s - l)!!(2s - 1)!! 
(n — s — l)\s\ 

immediately give Lemma ||. Thus the proof of Lemma || is complete. 

□ 

9. Proof of Theorem || 

If the side divided by the symmetry axis has odd length, the position s of the missing 
triangles ranges from 1 to n (see Figure ||). We can form the inner dual graph and denote it 
by G. Now we can proceed analogously to Section [3|, break the graph in two parts G + and 
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Lemmas QG, O and 12 



Figure 7. Af(H+(3,2)) equals M(i2+(2,2)). 




Figure 8. A tiling and lattice paths for R~ . 



G~ with the help of the matchings factorization theorem (see Lemma ||) . We convert G + and 
G~ back to regions R + and R~ of triangles and have to count rhombus tilings again. R+ and 
R~ are shown in an example in Figure |(| 
Thus, we have 

G = 2 n ~ 1 M(R+)M(R~). (11) 

Now we reduce the evaluation of M(R + ) to the evaluation of M(R + ), which we have already 

done in Lemma 0. M(R + ) and M(R + ) are related in the following way. The tiles of the upper 
half R + of the hexagon with sides n, n, 2m, n, n, 2m (as exemplified in Figure |]) sharing an 
edge with the border of length m are enforced as shown in Figure (the forced tiles are 
shaded). After removal of these tiles we are left with the upper half R + of the hexagon with 
sides n — 1, n — 1, 2m+l, n — 1, n — 1, 2m + 1. Thus, we have M(i? + (n, m)) = M(R + (n — l,m)) 
and Lemma ^ implies directly the following result: 



H(n + l) n 2 <»< J <n+i (2m + 2j 

n 5S (2j - 2) 



Af(12+) = — n+l J 1 ^. ( 12 ) 



The lower half R can be turned into a determinant in a manner analogous to Section ^ (see 
Figure g). 
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The starting and end points and the resulting determinant equal 
Ri - 



(2i — + m) for i ^ s 

(2s — 2, s + m — 1) for i = s 

Sj = (n + j- l,j - 1) 



if n + m — i 



A 



ij 



2 \m + i — j + 1 
n + to — s + 1 

m + s — j 



+ 



n + m — i 
m + i — j 



for j ^ s 
for i = s. 



Since the original problem and the claimed final result are invariant under n + 1 — s — > s, we 
can assume s ^ n. (The case n = s = 1 is trivial to check.) Then 



.4 



1 



m 



+ 



m 

to + n — j 



2 \m + n — j + 1 

{1 for j = n, 
else, 

since to is a nonnegative integer. It is easily seen that 

Aij(n, to, s) = Aij(n — l,m + 1, s — 1) for i, j < n, 

where Ay is defined in equation (|j). We expand deti<jj< n (^4jj(n, to, s)) with respect to row 
n and get 

det (Aij(n, m, s)) = det — 1, m + 1, s — 1)), 

Hence, 

M(R-(n,m,s)) = M(R7(n - l,m + l,s - 1)). 

Thus, Lemma ^ yields 

2( n 2 2 )- 1 H(n - l)(2n - 2s - l)!!(2s - 3)!! 



M(R~ 



(n- a )!(.-l)!lC (« + !)! 



n— 3 

n 

fc=i 



to + 1 + k + - 



min(fc,n-2-fc) n-1 



(to + 1 + k) 



min(fc+l,n— k) 



k=0 



Now we substitute the results of equations ( |l~2"l ) and (13) in (11). We get 

M(G) = 2 n - l M{R+)M(Rr) 

_ 2("^ 1 )" 1 H(n - 1) H(n + l)(2n - 2s - l)!!(2s - 3)!! 



(13) 



(14) 



(n-s)!( S -l)!n^ = o(2i)!nr=o + l 



n-2 



fc=2 



[] TO + /C + - 



min(fc— l,n— fc— 1) n 



Y[ ("I + fc) min ( fc ' n - fc+1 ) ]"f (2to + 2j 

2<i<j<n+l 



fc=l 



which can easily be transformed to the expression in Theorem 
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□ 



References 



[1] Mihai Ciucu, Enumeration of perfect matchings in graphs with reflective symmetry, J. Combin. Theory 
Ser. A 77 (1997), 67-97. 

[2] G. David and C. Tomei, The problem of the cahssons, Amer. Math. Monthly. 96 (1989), 429-431. 

[3] M. Fulmek and C. Krattenthaler The number of rhombus tilings of a symmetric hexagon which contain a 

fixed rhombus on the symmetry axis, II, preprint. 
[4] I.M. Gessel and X. Viennot, Determinant, paths and plane partitions, preprint, (1989). 
[5] C. Krattenthaler, An alternative evaluation of the Andrews-Burge determinant, in: Mathematical Essays 

in Honor of Gian-Carlo Rota, B. E. Sagan, R. P. Stanley, Progress in Math., vol. 161, Birkhauser 
[6] C. Krattenthaler, Generating functions for plane partitions of a given shape, Manuscripta Math. 69, 

(1990), 173-202. 

[7] P.A. MacMahon, Combinatory Analysis, vol. 2, Cambridge University Press, (1916); reprinted by Chelsea, 
New York, (1960). 

[8] J. Propp, Twenty open problems on enumeration of matchings, manuscript, (1996). 
[9] L. J. Slater, Generalized hypergeometric functions, Cambridge University Press, Cambridge, (1966). 
[10] J. R. Stembridge, Nonintersecting paths, pfaffians and plane partitions, Adv. in Math. 83 (1990) 96 — 131. 

Institut fur Mathematik der Universitat Wien, Strudlhofgasse 4, A-1090 Wien, Austria. 
E-mail address: teisenko@radon.mat .univie . ac . at 



16 



